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ABSTRACT 

Photoabsorption spectrum of the Xe@Cgo endohedral fullerene has been 
studied using the time-dependent-density-functional-theory (TDDFT), which 
represents the dynamical polarizability of an interacting electron system by 
an off-diagonal matrix element of the resolvent of the Liouvilliam superop- 
erator and solves the problem with the Lanczos algorithm. The method has 
been tested with the photoabsorption cross sections of the free Xe atom and 
Ceo fullerene. The result for the Xe@Ceo confirms the three main peaks 
observed in the recent measurement in the energy region of the Xe Ad gi- 
ant resonance and demonstrates the underestimation of the photoionization 
cross section by the experiment. It is suggested to use the current theoretical 
result around 22 eV to calibrate the experimental data. 
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1 Introduction 



Photoabsorption and photoionization cross sections of the Xe@Ceo molecule 
have been studied extensively by many groups both theoretically [H El El HI E] 
and experimentally [6J. Puska and Nieminen [JJ modeled for the first time 
the Ceo molecule using a spherical shell with an attractive potential and 
used a jellium shell to calculate the 4d photoabsorption cross section of a 
Xe atom inside the Cgo molecule. Their result demonstrated some oscilla- 
tions in the cross section curve. Amusia and co-workers [2j [3j S] used a 
5-like potential for the Ceo shell and the random-phase approximation with 
exchange (RPAE) method to describe the photoionization process. Their 
results indicate that the giant resonance is distorted substantially, exhibiting 
four prominent peaks [3] . Dolmatov and Manson [5J showed that two of those 
peaks disappear if the finite thickness of the Ceo cage is utilized. However, 
their result is not supported by the recent experiment [BJ. The experimen- 
tal data demonstrate strong enhancement of the photoionization cross sction 
as well as an oscillatory structure of the amplitude. The complicated ex- 
perimental resonance structure is also quite different from the calculation of 
the time-dependent density-functional theory (TDDFT) [7], which considers 
all the valence electrons to form a delocalized charged cloud and treats the 
residual ion-core as a classical jellium shell. Recently, Chen and Msezane 
performed a calculation on the confinement resonance of the Xe@Ceo- They 
solved the Schrodinger equation in three separate regions, r < r i: T{ < r < r Q 
and r > r Q using both regular and irregular solutions, where and r Q are, 
respectively the inner and outer radii of the Cgo spherical shell. The RPAE 
method was utilized to describe the photoionization process. This method 
has been successfully used in the study of the photoionization of the encap- 
sulated Ce 3+ [H [9] and Xe + [10] ions in Ceo- For the photoionization of 
the Xe@C6o molecule they obtained an improved agreement between theory 
and experiment [IT] . Although the photoionization of the Xe@Ceo molecule 



has been studied using various theoretical models, however, a first principle 
calculation is still needed. 

In this paper we have performed a time-dependent density-functional- 
theory (TDDFT) calculation for the photoabsorption spectrum of the Xe@Ceo- 
The calculations of the optical spectra are carried out in the linear response 
formalism. The most widespread calculation of the electronic excitation in 
TDDFT is the Casida's approach [12] . which reformulated the calculation of 
the response function into a generalized Hermitian eigenvalue problem. How- 
ever, this method is not suitable for calculating the spectra in a broad energy 
range. Recently a method has been developed, which determines the dipole 
susceptibility based on the application of the Lanczos algorithms to the first 
order expansion of the time-dependent Kohn-Sham equation in the linear re- 
sponse regime [T3J [HI [15] . Here the absorption spectrum is calculated from 
the imaginary part of the dipole susceptibility. The method avoids the ex- 
plicit calculation of the unoccupied states of the ground state and allows the 
entire spectrum of a molecule to be computed. This new implementation of 
the TDDFT has been used in the present Xe@Ceo calculation. 

2 Method and results 

We used the DM0I3 software of density-functional theory to determine the 
structures of the Ceo and Xe@C6o and the TDDFT to evaluate the photoab- 
sorption spectra. 

The TDDFT method was used first to calculate the photoabsorption spec- 
tra of the Ceo fullerene. The Cgo geometry optimization was performed us- 
ing the generalized gradient approximation (GGA) to the density-functional 
theory, with PBE exchange-correlation functional [TB] along with all elec- 
tron double numerical plus polarization (DNP) basis sets as implemented 
in the DM0I3 package [17J. The optimization of atomic positions proceeded 
until the change in energy was less than 5xl0 -4 eV and the forces were 



less than 0.01 eV/A. The resultant structure has the length of the shorter 
bond (hexagon- hexagon) 1.392 A and the longer bond (hexagon-pentagon) 
1.448 A. These value are in good agreement with the x-ray data value of 
1.391 A for the shorter bond and 1.455 A for the longer bond [18]. The 
structure was validated using a plane-wave approach as implemented in the 
QUANTUM-ESPRESSO [19]. A supercell of 18 A was constructed to elimi- 
nate the interactions among cells. An ultrasoft pseudopotential of the Rappe- 
Rabe-Kaxiras-Joannopoulos (RRKJ) [20] type and a kinetic energy cutoff of 
408 eV for the wave function and 2448 eV for the densities and potentials 
were employed in a standard ground state DFT calculation, yielding the 
Kohn-Sham eigenvalues and eigenvectors for all the occupied states. The 
dynamical polarizability of an interacting electron system is represented as 
an appropriate off-diagonal matrix element of the resolvent of a Liouvillian 
superoperator. The resolvent of the Liouvillian is evaluated through an al- 
gorithm based on the Lanczos method [T31 [LU E3- 

Fig 1 presents the photoabsorption spectrum for the Ceo after 3000 Lanc- 
zos iterations. Solid curve is from our TDDFT calculation and the dashed 
line reprsents the measurements [211 1221 [23]. The first three peaks at 3.48, 
4.36, and 5.36 eV agree very well with the other TDDFT and local orbitals 
calculation |24j. Their results show the peaks at 3.5, 4.4, and 5.4 eV. The 
Ceo fullerene has a giant resonance, which corresponds to a collective oscilla- 
tion of delocalized electrons relative to the ionic cage. The peak around 5-7 
eV is the 7r plasmon resonances. The largest peak around 22 eV is related 
to the a surface plasmon resonance, which was first predicted by the linear 
response theory [25]. The second collective resonance near 30-40 eV is as- 
sociated with a dipole-excited volume plasmon [20] ■ Our solid curve shows 
the second volume plasmon resonance predicted by the experiment [26] . The 
over all excellent agreement has been obtained between our TDDFT result 
and the measurements. 

The TDDFT method using the Liouville-Lanczos approach has also been 




Figure 1: (Color online) Photoabsorption spectrum of Cqq. Solid curve 
represents the result of the TDDFPT calculation. Dashed line are the ex- 
perimental data [2T1 1221 [23] 



used to evaluate the photoabsorption spectra for a free Xe atom. In order 
to study the 4d giant resonance of a Xe atom, an ultrasoft pseudopotential, 
which includes 4d electrons as valence electrons has been created using the 
ldl.x program of the QUANTRUM-ESPRESSO with the Xe ground state 
as the reference electron configuration. The pseudopotential has multiple- 
projector with PBE of GGA exchange-correlation functional and includes 
nonlinear core corrections. Fig. 2 shows the calculated photoabsorption 
spectra. Solid line is the calculated absolute differential oscillator strength of 
the photoabsorption of the Xenon atom in the energy region of 8.43 eV (5p 
threshold) -140 eV. Triangles are the experimental data obtained by Chan et 
al [27] using a low-resolution spectrometer. Solid squares and solid circles are 
respectively, the measured photoionization cross sections of the Ad subshell 
of Xenon atom by Holland et al [28] and Kammerling et al |29j. Note that 
the unit of photoionization cross section on the right vertical coordinate is 
the Mb, which is related to the differential oscillator strength through, 1.0975 
Mb=lxl0~ 2 differential oscillator strength unit (1/eV) [30]. Our calculation 
describes those measured data very well. 

The structure of the Xe@C 60 was obtained by adding a Xe atom into the 
center of the Cqo fullerene and then performing the same geometry optimiza- 
tion as that of the Ceo- The kinetic energy cutoff of 1020 eV for the wave 
function and 4080 eV for the density were used in the ground state calcu- 
lation. Fig. 3 shows the photoabsorption spectrum of the Xe@C 60 . In the 
energy range of the Xe Ad giant resonance, we can find several small wiggles 
which do not exist in the Cgo case. The enlarged details of these wiggles are 
presented in Fig. 4. Fig. 4 also shows the convergence situations of the 4d 
resonance. Dotted line, dashed line and solid line represent, respectively the 
1500, 3500 and 4500 Lanczos iterations. 

Fig. 5 is the comparison of our TDDFT calculation with the experiment 
in the energy range of the Xe Ad giant resonance. Solid circles with the error 
bars are from experiment [6]. Solid curve represents our calculation. Since 
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Figure 2: (Color online) Photoabsorption spectrum of free Xe atom. Solid 
curve is our TDDFT results. Triangles, solid squares and solid circles are 
respectively, from the measurements of [27], [28] and [29] . 
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Figure 3: Photoabsorption spectrum of Xe@C 60 
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Figure 4: (Color online) Comparison between the spectra at different iter- 
ations. Dotted, dashed and solid curves represent, respectively 1500, 3500 
and 4500 iterations. 
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Figure 5: (Color online) Comparison between the TDDFT calculation and 
the experimental data. Solid circles with error bars are taken from Ref. 
[6]. Solid curve is our TDDFT calculation. Experimental data have been 
enlarged seven times and shifted upwards by 0.25/eV. The calculated curve 
has been shifted to the right by 7 eV. 

the experiment underestimates the photoionization cross sections we have 
enlarged the measured cross sections by a factor of seven and shifted them 
upwards by 0.25/eV. The calculated Xe@C 60 curve is shifted to the right 
by 7 eV to match the three main peaks of the experiment. Our calculation 
confirms the three main peaks in the measurement and demonstrates that the 
experiment has underestimated the photoionization cross section of Xe@Cgo- 
We recommend the use of the TDDFT result at the low energy region to 
calibrate the measurement. For example, at 21.47 eV the photoionization cross 
section from TDDFT is 1633 Mb. 



3 Conclusion 



In conclusion, photoabsorption spectra for the Xe@Ceo molecule has been 
studied using the TDDFT method. The result confirms the three main peaks 
observed in the experiment [6J. The calculation also reveals the underestima- 
tion of the photoionization cross section by the measurement. It is suggested 
to use the TDDFT result around 22 eV to calibrate the measurement. 
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